Zurek-Kibble domain structures: The Dynamics of Spontaneous Vortex formation in 

Annular Josephson Tunnel Junctions 
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Phase transitions create a domain structure with defects, 
that has been argued by Zurek and Kibble to depend in a 
characteristic way on the quench rate. In this letter we present 
an experiment to measure the ZK scaling exponent a. Using 
long symmetric Josephson Tunnel Junctions, for which the 
predicted index is a = 0.25, we find a — 0.27 ±0.05. Further, 
we agree with the ZK prediction for the overall normalisation. 
PACS Numbers : 11.27. +d, 05.70.Fh, ll.10.Wx, 67.40.Vs 

Because phase transitions take place in a finite time, 
causality guarantees that correlation lengths remain fi- 
nite. Order parameter fields become frustrated, and de- 
fects arise so as to mediate the correlated regions with 
different ground states. Since defects are, in principle, 
observable, they provide an excellent experimental sig- 
nature for the way in which transitions are implemented. 

For condensed matter systems, which include the long 
annular Josephson Tunnel Junctions (JTJs) that we shall 
discuss below, Zurek jl]j|] suggested that causality alone 
is sufficient to determine the initial density of defects 
arising in a continuous transition. In this he paralleled 
proposals made by Kibble || in the context of quantum 
field theory models of the early universe. 

As applied to JTJs, the idea is as follows. Consider 
a thin linear JTJ with critical temperature T c , cooled 
through that temperature so that, if T(t) is the tem- 
perature at time t, then T(0) = T c . T (0) = —T c /tq 
defines the quench time tq. Suppose, at time t, that 
c(i) = c(T(t)) is the Swihart velocity vanishing at 

t — 0, and that £ a d(t) — Cad(T(t)) is the adiabatic heal- 
ing length (the Josephson length Xj(T(t))), diverging at 
t = 0. The first time that fluxons (or Josephson vortices), 
the defects of linear JTJs, can appear is at time i, when 

\L4t)\^c(t). _ 

For the case in hand, t has the form t = ^/tqtq, where 
r o T Q is the relaxation time of the longest wave- 
length modes. Details are given in papers ||,|tJ by two 



*E-mail: roberto@sa.infn.it 
^E-mail: myg@fysik.dtu.dk 

tE-mail: r.rivers@ic.ac.uk Permanent address: Blackett 
Laboratory, Imperial College, London, SW7 2BZ, U.K. 



of us (R.M and R.R). As a result, t q > t > r . If 
€a,d(t) ~ S,o{t/T Q y / for t ~ 0, where £ is simply re- 
lated to £,ad{T = 0), the cold fluxon size, then the initial 
domain size and fluxon separation is predicted to be 



£~£ad(*) =£0 



(1) 



where a = v/2. We term a the Zurek-Kibble (ZK) char- 
acteristic index. 

The arguments are not specific to JTJs. Prior to our 
experiment, five other condensed matter experiments had 
been performed to test the prediction (|l|) for the separa- 
tion of defects at their production, two experiments [|]f]] 
on superfluid 3 He, two |^p| , ^l| on superfluid 4 He, and 
one |l2| on high temperature superconductors (HTSC). 
In addition, an experiment |]l3| on JTJs by two of us 
(R.M and J.M) was compatible with ([!]), although it had 
not been performed with a test of (JlJ) in mind. 

Before discussing our new experiment a few comments 
are in order. The experiments ^|J|] on superfluid 3 He — B 
rely on the fact that, when it is bombarded with slow neu- 
trons, n+ 3 He — > p+ 3 H+760keV. The energy released in 
such a collision leads to a hot spot which, when cooled by 
its environment below T c , leaves behind a tangle of vor- 
tices (the topological defects in this system), tq is fixed 
by the nuclear process that breaks up the 3 He atom. 
With only a single data point conflating both normali- 
sation and a it is not possible to confirm the predicted 
value a — 1/4. However, both experiments are highly 
compatible with (Q), agreeing to a factor of a few in each 
case. 

In principle, the 4 He experiments 110], ll| , which use 



a pressure quench with a varying timescale tq to imple- 
ment the transition, allow for a more complete test. Yet 
again, vortices are the relevant defects. In practice, the 
most reliable experiment sees no vortices. This is 
not necessarily a sign of failure in that it has been sug- 
gested Q that the vortices decay too fast to be seen. 
This is irrespective 15 of whether a pressure quench, 
which preserves high thermal fluctuations, would of itself 
lead to somewhat different predictions. In this context, 
the vortices seen in an earlier 4 £fe experiment ]l0| ] were 
most likely an artefact of the experimental setup. 
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The fifth experiment [12], on high-T^ superconductors, 
measures total flux through a surface carried by the 
Abrikosov vortices. The vortex separation of ([!]) can be 
converted into a prediction for the flux, but no flux is seen 
in contradiction with this prediction, despite the phase 
separation that leads to the result being demonstrated 
elsewhere fTtj] . In this case there is no obvious explana- 
tion of the null result, despite later work jTjj that takes 
the effect of gauge fields into account fully. 



where J C (T) is the critical Josephson current at temper- 
ature T. Typically Aj(0) is in the 10 — 100 fj,m range and 
was equal to 7 fim for the sample discussed below. In 
Ref. M £o has been inferred as 



In summary, these early experiments (including |13|) 
have either provided one data point for ([!]), or have 
been null. Subsequently, two experiments of a very 
different type have been performed that permit vary- 
ing quench rates and so an estimate for a. The most 
recent | fl8| involves the Benard-Marangoni conduction- 
convection transition. The defects here are not associ- 
ated with the line zeroes of an order parameter field, 
and the viscosity-dependent a docs not match the ZK 
prediction, most likely for that reason. The more rele- 
vant experiment [ fL9| is carried out in a non-linear optical 
system, with complex beam-phase the order parameter. 
Increasing the light intensity (the control parameter in 
this case) leads to pattern formation (defects) at a crit- 
ical value. The predicted a = 1/4 is recovered to good 
accuracy, but agreement with normalisation is not stated. 

Our experiment, whose methodology and results we 
outline below, is also one in which, by varying tq, we 
can measure and compare a with its theoretical value, as 
well as confirming overall scale. 

In its essence, we quench a long annular JTJ through 
its critical temperature and count such fluxons or Joseph- 
son vortices as appear. In its idealised form the annular 
JTJ consists of two dimensionally identical annuli of su- 
perconductors of narrow width, superimposed upon one 
another, separated by an insulating barrier in the same 
plane. In practice, it is sufficient for the lower supercon- 
ductor to be a planar substrate upon which the other 
annulus sits. The effective theory for fluxons in JTJs 
H is the sine-Gordon model with respect to the field 
<t> — 4>\ — 02, the difference in the phases of the complex 
order parameter fields in the separate superconductors. 
The Josephson vortices are then the sine- Gordon kinks. 

The JTJs in our experiment are symmetric, by which 
is meant that the electrodes are made of identical super- 
conducting material with the same energy gaps and the 
same T c . This is confirmed by seeing that a) there is no 
logarithmic singularity in the voltage-current characteris- 
tic at finite voltages and b) the temperature dependence 
of the critical current is linear as T approaches T c . The 
outcome is that HQ v = 1/2. Therefore, at the time of 
their formation the separation of fluxons is expected to 
be given by (0) with a = 1/4. 

In terms of the parameters of the JTJs, the Josephson 
length at temperature T is f| 



Co 



2efi d s aJ c (0) ' 



U(T) = Xj(T) = 



2e^d a J c (T) ' 



The parameter 3 < a < 5 is given in terms of the su- 
perconductor gap energy A(T) and T c . As for r , it 
is given as r = £o/ c O; where c defines the behavior 
c(i) = cq (t / tq) 1 / 2 oi the Swihart velocity for the system 
near T = T c . If the thickness of the two superconducting 
electrodes differs, the effective thickness d s is the har- 
monic mean of the individual thicknesses ||. 

Our samples are high quality, 500 fim long, 3 \im wide, 
Nb/Al — Al ox /Nb JTJs fabricated on 0.5 mm thick silicon 
substrates using the trilayer technique (SNEAP) in which 
the junction is realized in the window opened in a SiO 
insulator layer. Details of the fabrication process can be 
found in Ref. (^0|. For all samples the high quality has 
been inferred by a measure of the I-V characteristic at 
T = 4.2 K . In fact, the subgap current I sg at 2mVwds 
small compared to the current rise AI g in the quasiparti- 
cle current at the gap voltage V g , typically AI g > 35I sg - 
The gap voltage was as large as V g = 2.76 mV and the 
maximum critical current I c was larger than 0.55A/ g for 
the overlap type junction. Furthermore, the application 
of a strong enough external magnetic field in the barrier 
plane completely suppressed any Josephson structures in- 
dicating the absence of electrical shorts in the barrier. 
At an order of magnitude level, we have To ~ 0.1 ps and 
£o ~ 1 M m - As a result £ ~ 1 mm for tq ~ 1 s. We will 
be more specific later. 

Our JTJs have a critical temperature T c = 8.95 if, 
whereas the individual superconductors have a critical 
temperature of 9.1 K. Even at our fastest quench the 
conductors are superconducting, by which is meant that 
the Cooper-pair order parameter field has achieved its 
final magnitude, 1 ms before the JTJ can develop fluxons. 
This is necessary for (0) to be valid without modification 
since only then is <f> the relevant order parameter. 
In order to vary the quenching time over the widest 
possible range, we have realized the experimental setup 
shown schematically in Fig.l. The annular JTJ samples 
are fabricated on a chip (shown endways on), mounted 
to a Cu block by a thermally insulating teflon sheet. The 
entire system is enclosed in a vacuum-tight can immersed 
in the liquid He bath. In all cases the heat is removed 
from the system by He exchange gas, using a manual 
pump. By varying the pressure of the gas we can mod- 
ify the rate of cooling of the sample. On the other side 
of the block is a 50 f2 carbon resistor, that enables us 
to heat and cool the JTJs on a relatively long timescale 
(tq from about 1 to 10 s) depending on the He exchange 
gas pressure inside the can. With a — 1/4, we need to 
vary tq by at least two orders of magnitude. This is 
not possible just by heating the block alone (even if it 
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were smaller). To extend the range of tq we mounted 
a small, pulse-driven, surface 100 f2 resistor on the same 
side. This permits much smaller thermal cycles (0.07 to 
0.2 s). These two completely different quenching tech- 
niques provide timescale ranges that do not overlap, leav- 
ing a gap between 0.2 and Is, that would require a third 
quenching technique to be filled. 

The whole system is then enclosed in a /^-metal 
shielded cryostat. The temperature of the JTJ is mon- 
itored by measuring the junction gap voltage, which is 
proportional to the known superconductor gap energy 
A(T(f)). All quenches were taken from 10 A" to 4.2 K 
through T c = 8.95 if. By making use of the Thouless 
equation pjj] , it is possible to infer the temperature of 
the JTJ from the gap voltage in the range 8.2 K to 4.2 K, 
extremely accurately at the upper end of the range, and 
with fluctuations of a few per cent at the lower. Whether 
for slow or fast cooling an excellent fit to the temper- 
ature of the JTJ in this range, for initial temperature 
Ti n = 10 K and final temperature Tfi n = 4.2 K, is given 
by the thermal relaxation equation 



T(t) = T ln + (T in ~ Tfi n ) e" 



-(t-to)/r 



(2) 



where r is the relaxation time which sets the cooling 
timescale. See Fig. 2 for an example. This equation is 
then used to extrapolate dT(t)/dt to the vicinity of Tc 
and yields tq ~ 1.7 t. 

On cooling the system in this way, we expect fluxons 
to appear from the inhomogeneity of <p at the transition, 
according to (Q). In the absence of any current through 
the barrier or applied external magnetic field the fluxons 
are in indifferent equilibrium as far as the barrier is ho- 
mogeneous and pin-hole free. (In reality, there is a small 
pinning effect so that, after a short transient, the flux- 
ons are static.) To make them visible, and countable, 
we apply a bias current, whereby they move as mag- 
netic dipoles under the resulting Lorentz force, at speed 
less than the Swihart velocity. According to Josephson 
theory N travelling fluxons (and antifluxons) with speed 
v develop a voltage V = N<& v/C across the junction, 
where $o = h/2e is the flux in a (Josephson) fluxon 
and C is the annulus circumference. This voltage can 
be measured, and the fluxon number determined. A de- 
tailed description of the experiment and the data will be 
given elsewhere j22|. 

Many samples have been measured, but only one had 
such a large critical current density (and sufficiently small 
Xj(T)) that only 3000 thermal cycles were enough to 
get reliable statistics, and it is this sample that we shall 
discuss now. However, we stress that, within the less 
good statistics of the other samples, none gave results 
that were incompatible with (Q). 

The symmetric annular JTJ with which the experi- 
ment was performed had a circumference C — 500 /im, 
and width Ar = 4 /im. The effective superconductor 
thickness was d s rs 250 nm. At the final temperature 
Tfin = 4.2 AT, the critical current density was J c {Tfi n ) — 



3050 A/ cm 2 , the quality factor was V m = 49mT^and the 
Josephson length was \j(Tf in ) — 6.9 /im. The velocity 
co is calculated to be Co ~ 2.2 x 10 7 m/s. From this, we 
infer that £o ~ 3.8 /im and tq w 0.17ps. It follows that 
£ ~ 5.9mm for tq = Is, several times the circumference 
of the JTJ. 

As a result, the likelihood of finding a single fluxon is 
small. We estimate the probability of finding a fluxon in 
a single quench to be 



(3) 



where, from (0), cr = 0.25. 

In Fig. 3 we show a log-log plot for the measured Pi as a 
function of tq , changed by varying the exchange gas pres- 
sure and by using both the fast and slow quenching tech- 
niques. Each data point corresponds to 300 quenches. 
Although the data does not distinguish between a single 
fluxon and a fluxon plus an antifluxon pair, the likeli- 
hood of the latter is sufficiently small that it can be ig- 
nored. Similarly, the data does not distinguish between 
no fluxons and a fluxon-antifluxon pair, with similar er- 
rors. We observe that the points are quite scattered, 
meaning that the data are statistically poor. Further, 
for the reasons given earlier there is a gap between fast 
and slow quenches, with the possibility for systematic 
error. Nonetheless, we have clear evidence that i) the 
trapping of a fluxon occurs on a purely statistical basis 
for identical conditions of each thermal cycle and ii) the 
probability to trap one fluxon is larger when the tran- 
sition is performed at a faster speed (smaller quenching 
time) in accordance with the causality principle. This 
complements our qualitative results from other samples 
with smaller J c (0) (for which the statistics is too poor to 
permit a fit to O)) that, at fixed tq, the probability of 
finding a fluxon decreases with increasing £o- 

Regardless of the data spread, to test (pf) we attempted 

-b 



to fit the data with an allometric function Pi 



Q . 



with a and b being free fitting parameters. We found, 
for the coefficient a, the best fitting value of 0.1 ± 10% 
(taking tq in seconds). This is in excellent agreement 

with the predicted value of Ctq^ 4 /£o ~ 0.08s 1 / 4 , given 
the fact that we only expect agreement in overall normal- 
ization to somewhat better than an order of magnitude. 
After the failure of the experiments discussed in (h] [j^] 
to find (reliable) defects at expected densities, if at all, 
this experiment shows that the ZK estimate remains sen- 
sible. Further, the best fitting curve, shown by the solid 
line in Fig. 3, has a slope b = 0.27 ± 0.05, in remarkable 
agreement with its predicted value of 0.25. 

The ZK scenario needs further testing with JTJs for 
which there is a greater likelihood of observing fluxons. 
In H and H we observed that this is the case for sig- 
nificantly non-symmetric JTJs, for which the value of a 
is a = 1/7; a further experiment, with markedly non- 
symmetric JTJs, is being planned. 
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I. FIGURE CAPTIONS 

Fig.l Sketch (dimensions are not to scale) of the cryo- 
genic insert developed to perform the junction ther- 
mal cycles with a time scale changing over a broad 
range. The JTJs are fabricated within the chip 



(shown endways) , to which the surface mounted re- 
sistor (SMR) is attached. The whole is surrounded 
with liquid He. 

Fig. 2 Time dependence of the junction temperature dur- 
ing a "slow" thermal cycle from 10 K to 4.2 K, as 
inferred from the junction gap voltage. Only be- 
low the horizontal dashed line does the continuous 
curve describe the temperature. The thick dashed 
line is the best fitting curve to (||), to be extrapo- 
lated to T c . 

Fig. 3 Log-log plot of the measured probability Pi to trap 
one fluxon versus the quenching time tq. Each 
point corresponds to 300 thermal cycles. The solid 
line (slope b = 0.27) is the best linear fit, in good 
agreement with the 0.25 value expected for sym- 
metric JTJs. Errors in tq are a few per cent, and 
systematic errors in Pi due to the neglect of fluxon- 
antifluxon pairs are O(Pl), again a few percent. 
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